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Microbubbles, either in the form of free gas bubbles surrounded by a fluid or encapsulated bubbles used 
currently as contrast agents for medical echography, exhibit complex dynamics under specific acoustic exci¬ 
tations. Nonetheless, considering their micron size and the complexity of their interaction phenomenon with 
ultrasound waves, expensive and complex experiments and/or simulations are required for their analysis. The 
behavior of a microbubble along its equator can be linked to a system of coupled oscillators. In this study, the 
oscillatory behavior of a microbubble has been investigated through an acousto-mechanical analogy based on a 
ring-shaped chain of coupled pendula. Observation of parametric vibration modes of the pendula ring excited 
at frequencies between 1 and 5 Hz is presented. Simulations have been carried out and show mode mixing 
phenomena. The relevance of the analogy between a microbubble and the macroscopic acousto-mechanical 
setup is discussed and suggested as an alternative way to investigate the complexity of microbubble dynamics. 

PACS numbers: 43.25.Ts,43.20.Wd,43.25.Gf 


I. INTRODUCTION 

When subjected to an external acoustic field, bubbles 
can undergo complex radial oscillations. This oscilla¬ 
tory behavior has been an important and continuously 
developing research subject since the beginning of the 
twentieth century. The investigation of bubble dynamics 
started with the work of Lord RayleiglP^, who was man¬ 
dated by the Royal Navy to explain the origin of dam¬ 
ages on submarine propellers. Rayleigh focused on the 
oscillatory behavior of cavitation bubbles suspended in a 
fluid. He showed that the overpressure generated by the 
oscillations and the collapse of bubbles could explain the 
damages caused on propellers. In the 30’s, Minnaert was 
interested in the origin of the sound of running watei!^. 
He supposed that bubbles oscillating periodically in wa¬ 
ter were at the origin of rivers whispering. After these 
seminal works, many studies on bubble oscillations have 
been carried out in different research fields. 

In addition to the radial motion, bubbles can show non- 
spherical oscillations or vibration modes. These vibration 
modes, characterized by an index n, were first analyzed 
theoretically at the interface between immiscible and in¬ 
compressible fluids with spherical symmetric^. Later, 
NeppiraJ^ analyzed the acoustic response from gas bub¬ 
bles suspended in a fluid and subjected to sound fields. 
Eller and CrunP, and Prosperetti et alW focused on the 
instability of the motion and the nonlinear dynamics of 
a bubble within a sound field. In the 90’s, the discov¬ 
ery of sinde bubble sonoluminescence (SBSL) by Gaitan 
and GrunP, led to additional studies on nonlinear oscilla¬ 
tions of bubbleFinally, with the use of ultrasound 


contrast agent the understanding of bubble dynamics 
have found a renewed interest in t he fie ld of ultrasound 
imaging and targeted drug deliver}EE^. Therefore, the 
problem of vibration modes in bubbles is still under in¬ 
vestigation. 

In this study, the interaction between ultrasound and a 
microbubble, and especially the appearance of vibration 
modes, are studied using a macroscopic analog system. 
The mechanisms underlying its nonlinear behavior are 
sought to improve our understanding of the microbubble 
dynamics. The study of a single microbubble is a dif¬ 
ficult task, particularly because of the interaction with 
other microbubbles or microstreaming that sweep away 
the bubble. Moreover the smallness and the complexity 
of the phenomena involved, require complex modeling 
and expensive experiments. In this work, we propose the 
use of a macroscopic mechanical analog as an alternative 
way to investigate microbubble dynamics. 

The study of analog models in the field of physics is 
a tool that allows recreating in the laboratory phenom¬ 
ena that are difficult to observe directly. Goncerning 
microbubbles, some of its dynamic features (vibration 
modes, subharmonic oscillations, chaos) can be captured 
by a system of coupled oscillators. Here, the oscillatory 
behavior of a microbubble and its vibration modes are 
investigated through the use of a macroscopic analogy 
consisting of a chain of coupled pendula, parametrically 
excited by a vertical force. Based on a discrete nonlinear 
model of coupled pendula and its continuous limit de¬ 
scribing low frequency excitations, vibration modes are 
investigated theoretically and experimentally. This ap¬ 
proach is used here to set up a formal basis for the 
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acousto-mechanical analogy of a gas microbubble. 

The structure of the paper is as follows: Mathematical 
formulations of a gas microbubble in an ultrasound field, 
and the analog acousto-mechanical system are given and 
compared in sections II and III respectively. The ex¬ 
perimental setup and the results of the measurement of 
vibration modes in the macroscopic system are described 
in section IV, where numerical simulations are also pre¬ 
sented. The conclusions of the study are presented in 
section V. 


II. MICROBUBBLE DYNAMICS 


In the presence of an acoustic field, microbubbles can 
be forced to oscillate in different ways. The most common 
oscillation mode is the radial mode (Fig. [^, where the 
bubble compresses and expands radially, maintaining its 
spherical shape. The basic model describing the radial 
dynamics of a bubble is the Rayleigh-Plesset equation for 
the time dependent radius R{t) 

p (^Rii + =Pg-Po- PA{t) (1) 

with p the density of the surrounding fluid, Pg = 
(^0 ^ pressure inside the bubble, Pq 

is the hydrostatic pressure, Rq is the equilibrium radius 
of the bubble. Pa the acoustic pressure, a the surface 
tension of the bubble, p the dynamic viscosity and 7 the 
polytropic exponent. Some generalizations of this model 
have been proposecP^HIll^ The resonance frequency of 
the radial mode of the bubble, the so called Minnaert 
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FIG. 1. A bubble oscillating volumetrically within an acoustic 
field. For positive and negative pressures, compression and 
expansion phases are observed respectively. 


frequency, can be obtained from Eq. 0 and is given by 


/ 


27rRo VP/ 


( 2 ) 


Under typical conditions, the relation /Rq is a mag¬ 
nitude of order 1. In particular, for a bubble in water 
at standard pressure (Po=100 kPa, p=1000 kg/m^), this 
equation gives the condition /Rq ^ 3.26 m/s. 

Microbubbles can also undergo non-spherical oscilla¬ 
tions (see Fig.[^ left column) through instabilities at the 
gas/fluid interface. In this case, the radius becomes a 
space-dependent function. 


R{t) ^ R{9,<p,t) = R{t) + (3) 

where R{t) describes the evolution of the radial mode 
and ^{O^cpR) a perturbation depending on the spherical 
coordinates 0 and 0. As Fig. (left) shows, microbub¬ 
bles can develop different surface patterns that depend 
on the excitation parameters (amplitude and frequency 
of the ultrasound wave) and the bubble radius. The ra¬ 
dial oscillation corresponds to the mode n = 0; the mode 
n = 1 corresponds to the displacement of the center of 
mass. Non-spherical surface modes are those modes with 
n > 2. 

To describe non-spherical modes in bubbles, a com¬ 
mon analytical approach is to expand the perturbation 
of the radial mode, (^(P,(/9, t), on the basis of spherical 
harmonic^ 


ae,‘P,t) = J2an{t)Y:^{e,A, ( 4 ) 

n,m 


where an{t) is the time-dependent amplitude of the sur¬ 
face mode of index n, and Y^{0^ ip) is the spherical har¬ 
monic defined as 


( 5 ) 

where P^ are Legendre polynomials. Although (p) 

defines a large set of possible surface modal oscillations, 
experiments^ show that the observed modes present 
symmetry along the axis of the incident ultrasound beam, 
corresponding to m = 0, also known as zonal harmonics. 
Then, spherical harmonics V^(6>,0) reduce to Legendre 
polynomials Pn{cos0). The evolution equation for the 
amplitude of each mode can be found by matching veloc¬ 
ity potentials and pressures at both sides of the interface, 
and linearizing for small amplitude^. 


3R 2(n + 2)(2n + l) 


R ^ 


pR^ 


p \ an Y 


(n + l)(n + 2)(7 2{n + 2)pR R\ 
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FIG. 2. Representation of microbubble surface modes with n 
the mode’s order. The left pictures are experimental results 
from Versluis et a/.^, and show a selection of surfaces modes 
observed for different microbubble radius between 36-45 /xm. 
The two columns in the center are 3D analytical solutions 
from Eq. and its equatorial cross-section. The right col¬ 
umn shows the corresponding vibration modes experimentally 
observed with the acousto-mechanical system, as described in 
Section IV. 


where, /x is the viscosity and a the surface tension. Defin¬ 
ing h{t) = the equation can be simplified a^ 


/ (n — l)(n + l)(n + 2)cr 3R‘^ {2n-\-l)R\ 

4^ 2^ J = 0 


( 7 ) 

Within this approach, each mode n obeys to the equa¬ 
tion of a harmonic oscillator, with time dependent coef¬ 
ficients. The resonance frequencies of the surface modes 
readily follow from Eq. 0 . by considering the static con¬ 
dition R = Rq, and R = R = 0. 


CUn 



l)(n + l)(n + 2 )c7 
pRl 


(8) 


which is the Lamb expression for surface modes for a free 
gas bubble. 

The acoustic pressure term Pa in Eq. 0 is usu¬ 
ally a harmonic function with angular frequency cJe, i.e. 


= Pa cos(co’et). For sufficiently small amplitudes 
Pa, the bubble response will be also harmonic at the 
same frequency, i.e. 

R{t) = Rq p Rs cos(cjet), (9) 

with <C i?o- Substituting Eq. 0 into Eq. 0 and lin¬ 
earizing, a Mathieu equation for each surface mode is 
obtainecP^ 


bn + ^ cos(wet)^ bn = 0. 

( 10 ) 

The Mathieu equation is a special case of a lin¬ 
ear second-order homogeneous differential equation with 
time dependent coefficients, and appears in many appli¬ 
cations in physics and engineering, specially in the de¬ 
scription of parametrically driven system^, as the sys¬ 
tems considered in this work. 


III. THE PARAMETRICALLY DRIVEN CHAIN OF 
COUPLED PENDULA 


In this section the model equations for the macroscopic 
analog of the microbubble are formulated. First, we con¬ 
sider the exact problem of the discrete lattice of cou¬ 
pled masses, that corresponds to our experimental sys¬ 
tem. Later, the continuum limit of this model is used to 
establish the analogy with gas microbubble, by deriving 
an equation isomorphic to Eq. (10). 


The discrete lattice 


The equation of motion of a pendulum of length L 
given by 


IS 


0 + cjq sin© = 0, 


(11) 


where 0(t) is the angle with respect to the vertical, and 
ujQ = gjL^ g being the acceleration due to gravity. 
Here, for the sake of simplicity we have ignored dissipa¬ 
tion, which can be included in Eq. (11) by an additional 


term j3Q. When a set of oscillators are coupled to their 
nearest neighbors, they form a lattice or chain, support¬ 
ing waves. Consider that the lattice is subjected to a 
parametric forcing with displacement amplitude he and 
angular frequency cje, then its motion is described by 

©i + (^0 +(^coscJet) sinQi - c^(0i+i - 20^ + ©i-i) = 0, 

( 12 ) 

where 0^ is the angle of the i-th pendulum, g = AujlhelL 
is the forcing parameter, and c is a constant denoting the 
strength of the coupling, equivalent to the speed of sound 
for the waves for the lattice waves. Lattice models like 
Eq. (12) have been extensively studied, in the context 
of Frenkel-Kontorova chaiiP^l. Most of the work is theo¬ 
retical, and focused on the formation of localized states 
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can be obtained, for negligible damping and forcing, af¬ 
ter linearization and assuming a solution in the form of 
a harmonic wave, exj[)i{kna — ujt). This results in the 
relation 



(14) 


As shown in Fig. ib) , which illustrates the chain 
viewed from the top, the ring of oscillators is equivalent 
to an equatorial section of a bubble with equilibrium ra¬ 
dius i?o, 


Ri{t) = Rq LsmQi{t) (15) 

where Rq is the radius of the ring (equivalent to the equi¬ 
librium radius of the bubble). Assuming small angles 
sin0i 0^, the position of each pendulum at a time t 
can be given by Ri{t) = Rq LQi{t) 


B. Continuous description 


FIG. 3. Schematic representation of the chain of pendula 
coupled by V-shaped strings and knots, hanging from a rigid 
support, (a) Lateral view of a section of the chain, (b) top 
view of the whole chain in a circular arrangement. Deviation 
of one mass with respect to its equilibrium position, and its 
coordinates are shown for illustration. 


(breathers, solitons and kinks). There are few exper¬ 
imental work^^^H^, where different parametric driving 
mechanisms and coupling have been implemented. In 
Ref.l^, coupling is achieved by strings and knots, while 
in Ref.l^ torsional springs are used. In all the cases, a 
linear arrangement of pendula was considered. Boundary 
conditions typical of this configuration are either fixed or 
free ends. 

Here, in order to mimic a bubble-like behavior, we con¬ 
sider the chain in a circular arrangement. This implies 
the periodic boundary condition 07V+1 = 0i. In addition 
to this, the circular configuration incorporates curvature 
effects, which are not present in the linear chain. We 
assume an V-shaped coupling (see Fig. [^a)), for which 
coupling strength is given by 


The bubble surface is a continuum, while the chain of 
pendula is a discrete system. Therefore, analogies must 
be searched in the limit where a continuum version of 
Eq. (12) applies. This is the case when we restrict our 
analysis to long-wavelength modes (low index n), where 
the mode scale are much larger than the distance between 
two pendula, knCi « 1, or na/Ro « 1. Considering 
this limit, the discrete angle coordinates &i{t) can be 
replaced by the continuous function 0(x,t). 

The equation of motion Eq. ( [T^ leads to the paramet¬ 
rically driven, small amplitude, sine-Gordon model, that 
reads 


0 - C^Qxx + (^0 + ^ cos(cjet))0 = 0 (16) 


where 0 is a surface deformation related to the bubble 
radius as 0 = (R — Rq)/L. 

Solutions of Eq. (16) can be expressed as a superpo¬ 


sition of normal modes. Eor a linear (straight) chain, a 
proper basis is given by the harmonic functions 


<d{x,t) = y^^bnjt) cos{knx). 


(17) 


gda^ 

AL{L - d) ’ 


(13) 


where d is the distance from the knot to the support, 
and a the distance between masseP^. Note that, in this 
model, the coupling strength can be varied by selecting 
the position of the knot. 

The chain supports different oscillation modes, labeled 
with an integer index n, corresponding to standing waves 
with wavenumber kn . Eor the circular chain of radius Rq , 
the relation kn = ti/Rq holds. The relation between driv¬ 
ing frequency and mode index (the dispersion relation) 


Some solutions of Eq. have been discussed for a 

one-dimensional straight geometr^EU. 

In the ring geometry considered here, angular coor¬ 
dinates are more suitable to describe the position of a 
point on the deformed ring, and we use the transforma¬ 
tion X = RO, where 0 = [0, 27r[ and R = Rq is the equi¬ 
librium radius. The Laplacian operator takes the form 
^xx = R~^Oee- A proper basis for the expansion in this 
case is 


©(6», i) = E bn{t)Pn{cos 9), (18) 

n 
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where P^(cos6>) is a Legendre polynomial and 0 is the 
angular coordinate. 

The equation for the temporal evolution of the mode 
amplitudes can be obtained by substituting Eq. 0 or 
Eq. into Eq. and projecting over the differ¬ 

ent modes (Galerkin projection) using the orthogonality 
properties of harmonic functions or Legendre polynomi¬ 
als. Independently of the basis chosen for the expansion, 
the following equation is obtained 

K + + »? COs(Wet)) bn = 0, (19) 


which is a Mathieu equation, where the parametric exci¬ 
tation amplitude is and the frequency on the 

n-th mode is given by 


Z Z I z 

=^^0 +C 

Uq 


( 20 ) 


The same result is obtained from the exact dispersion 
relation Eq. (14) for the discrete system, evaluated in 


the limit of low index n, where na/Ro << 1. 


C. Analogy 


The previous theoretical study stablish a number of 
analogies between the ultrasound driven gas bubble 
and the ring of coupled pendula subjected to a time- 
dependent acceleration. Both systems are described by 
the same Mathieu equation, with different coefficients. A 
similar behavior is therefore expected for both systems, 
and similarity parameters are listed in table |lj Eor typical 
experiments, microbubbles {Rq 10“^ m) are insonified 
in the MHz frequency range. Note that, from Eq. (|^, the 
product fRo is fixed in microbubbles, but can be tuned 
in the mechanical macrobubble, where / and Rq can be 
varied independently. The pendula ring used here, with 
Rq ^ 0.3 m, is designed to be excited in the Hz frequency 
range. We can therefore consider / • Pq as a geometrical 
invariant. 


TABLE 1. Table of analogies between parameters of mi¬ 
crobubble and those of the pendula ring. 


Parameter 

Microbubble 

Pendula ring 

Ro [m] 

~ 10"® 

^ 0.3 

fe [Hz] 

~ 10® 

^ 3 

fe-Ro 

1 

1 


(n-l)(n + l)(n + 2)^ 


fv 

(211+1), .1 (rA 

2 

4-? (f) 


IV. EXPERIMENTS 
A. Experimental setup 

The setup consists of an aluminum ring, on which pen¬ 
dula of mass m = 6 g are fixed with nylon strings forming 
a “V” shape with the vertical axis as shown in Eig. [^a). 
The coupling between pendula is obtained by overlapping 
the strings and fixing them by a knofP^. The parameter 
with a stronger infiuence on the dynamics of the chain 
is the coupling strength. Here, we focus on a ring with 
a medium coupling. The ring has a radius P = 31 cm, 
= 54 pendula, pendulum’s length is P = 10 cm and 
the Idistance between the ring and the node is d = 5 cm. 
As Eig.|^ shows, the pendula ring lies on the excitation 
system. The mechanical excitation system consists of a 
subwoofer loudspeaker driven by an arbitrary waveform 
generator (Agilent33220A) through an audio amplifier. 
Thus, the pendula ring is attached to the subwoofer cone 
and therefore is excited mechanically by a vertical oscil¬ 
latory force, as shown in Eig.|^ The sinusoidal excitation 
varies from /g = 0 to 5 Hz with amplitudes varying from 
A = 0.5 to 3.5 Vpp, corresponding to a vertical displace¬ 
ments ranging from he = 0.5 mm to 3.5 mm. The pen¬ 
dula are driven near the volumetric resonant frequency 
and its double to allow the development of surface modes 
during the mechanical excitation through parametric in¬ 
stability. 

The motion of pendula has been recorded from the top 
with a video camera. The data processing has been per¬ 
formed with the software Image J and the plug-in MJ 
Track. This plug-in enable to track the motion of each 
pendulum and thus to determine the distance Ri{t) be¬ 
tween the center of the ring and a pendulum at a fixed 


camera 

□ 


pendula ring 



FIG. 4. Scheme of the experimental setup. The ring of cou¬ 
pled pendula is driven vertically by an sinusoidal force gener¬ 
ated by a loudspeaker. 
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FIG. 5. Experimental mode n = 3 (symbols) and the corre¬ 
sponding fit to a Legendre polynomial PsicosO). Two differ¬ 
ent representations, in cartesian (left) and polar (right) coor¬ 
dinates 


time t. Applying a spatial fast Fourier transform on 
Ri{t)^ the amplitude of each vibration mode of order n 
can be obtained. 


B. Results 


A set of measurements has been carried out for exci¬ 
tation amplitudes between he = 0.5 and 3.5 mm and fre¬ 
quencies varying from /e = 0 to 5 Hz. Vibration modes 
have been observed up to n = 22, including an unstable 
volumetric mode for n = 0. The mode n = 1 correspond¬ 
ing to the displacement of the center of mass has also 
been observed. An example of the observed patterns is 
given in Fig. for an excitation amplitude he = 2 mm 
and frequency fe = 3.20 Hz. Here, the pendula ring 
shows a mode n = 3. A fit to the corresponding Legen¬ 
dre polynomial, shown in continuous line, shows a good 
agreement. 

A phase diagram showing the domain of existence of 
vibration modes n for different excitation parameters is 
presented in FigJ^ The numerical simulation of Eq. (12) 
evidences the existence of different sets of resonances. 
Due to the damping, there exist a threshold for the ex¬ 
citation of the different modes. The lowest threshold is 
obtained for the so-called 2:1 or subharmonic resonance, 
where the frequency of the excitation is twice the natu¬ 
ral frequency of the pendulum. Whithin each resonance 
set, each mode has its own instability region, or Arnold 
tongues, which are represented in Fig. (|6) with a differ¬ 
ent color for each mode. The dark shaded region cor¬ 
responds to numerically unstable solutions. Note also 
the existence of a second set at lower frequencies, with a 
higher threshold, denoting the 1:1 resonance, that corre¬ 
sponds to modes excited when driving the s ystem at the 
natural frequency of the pendulum /o = {g/L) ^1.6 

Hz. The symbols in Fig.(|^ correspond to experimen¬ 
tal data. From red to blue, the mode number increases, 
similarly to the numerical results. No vibration modes 
are observed at frequencies below fe/fo = 0.95, as there 
are no unstability regions predicted by the theory. From 
fe/fo = 0.95 to 1.90, a first group of vibration modes 




fe (Hz) 


FIG. 6. (Top) Phase diagram showing the appearance of vi¬ 
bration modes for different excitation parameters. The sym¬ 
bols correspond to the experimental data whereas the colored 
areas are the instability regions obtained numerically. (Bot¬ 
tom) Oscillation frequencies of the pendulum (denoted by /o) 
versus the excitation frequency (/e). On both graphics, each 
symbol corresponds to a vibration mode. 


appears when the excitation amplitude is greater than 
a threshold. These modes oscillate at a frequency equal 
to the excitation frequency (/osc = /e), and correspond 
to the 1:1 parametric resonance of the system. From 
fe/fo = 1-90 to fe/fo = 3.17, a second group of vibra¬ 
tion modes oscillating at subharmonics of the excitation 
frequency (2 /osc = fe) is observed, corresponding to the 
2:1 resonance. Therefore, two sets of resonance modes 
have been observed. Comparing the trend of simulation 
and experimental data shown in Fig. we can conclude 
that there exists a good agreement. 

It is important to note that pure modes exist only 
in a narrow region. Close to the parametric instabil¬ 
ity threshold, the Arnold tongues corresponding to each 
mode n overlap. Thus, most of the modes observed actu¬ 
ally correspond to a mixing of neighboring modes, with 
one mode being clearly dominant. Mode mixing is there¬ 
fore expected to occur too in ultrasound-driven bubbles; 
however, such mode mixing has been postulated but not 
analyzed or observed in such microbubbles. The mechan- 
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ical analog presented here could be used to explore these 
complex phenomena. 

Finally, we note that the maximum excitation has been 
measured at /e = 5 Hz (///o = 3.2). Above this fre¬ 
quency, measurements cannot be performed due to the 
setup vulnerability. Oscillations becomes jerky and the 
coupling between pendula often breaks. 


V. CONCLUSIONS 

The interaction between a microbubble and an acous¬ 
tical field has been studied through the use of a macro¬ 
scopic acousto-mechanical analogy. As for real bubbles, 
vibration modes and modes mixing have been observed. 
It is known that, in the case of microbubbles, vibration 
modes display a strong subharmonic behavior. Here, 
with the pendula ring, vibration modes have also been 
excited in the region where /osc = /e, corresponding to a 
1:1 resonance. 

Localized modes like breather have not been discussed 
here but a previous study on the pendula ring shows 
the appearance of such rnode^^^^. One can thus ex¬ 
pect to observe these oscillatory behaviors in microbub¬ 
bles. These results allow to offer new insights in the 
study of microbubble’s dynamics. The analogy between 
the macroscopic behavior of the pendula ring and the 
microscopic behavior of a microbubble presented here is 
valid for low order modes (i.e. in the limit na/Ro 1) 
and provides important insights for the comprehension 
of microbubble dynamics. Further studies are under¬ 
way, particularly an expansion of the analogy to consider 
encapsulated microbubbles to get information for both 
imaging and therapeutic applications using contrast mi- 
crobubbles. 


ACKNOWLEDGMENTS 

Thanks to Jean-Marc Gregoire and Jean-Yves Tartu 
for their fruitful help in the design of the experimental 
setup. Jennifer Chaline was the recipient of a PhD fellow¬ 
ship from the Region Centre (France). Financial support 
from the University Francois-Rabelais doctoral school of 


Tours was granted to finance a research stay on this topic 
at the Universidad Politecnica de Valencia (Spain). A. 
Mehrem is recipient of a research fellowship from Gener- 
alitat Valenciana (Santiago Grisolia program). 

^L. Rayleigh, Philosophical Magazine 34, 94 (1917). 

^M. Minnaert, Philosophical Magazine 16, 235 (1933). 

^M. Plesset, J. Appl. Phys. 25, 96 (1954). 

Neppiras, J. Acoust. Soc. Am. 46, 587 (1969). 

^A. 1. Eller and L. A. Crum, J. Acoust. Soc. Am. 47, 762 (1970). 
®A. Prosperetti, L. A. Crum, and K. W. Commander, J. Acoust. 
Soc. Am. 83, 502 (1988). 

^D. F. Gaitan, L. A. Crum, C. C. Church, and R. A. Roy, J. 
Acoust. Soc. Am. 91, 3166 (1992). 

®M. P. Brenner, D. Lohse, and T. Dupont, Phys. Rev. Lett. 75, 
954 (1995). 

®F. Gaitan and G. Holt, J. Acoust. Soc. Am. 103, 3046 (1998). 
H. Augsddrfer, A. K. Evans, and D. P. Oxley, Phys. Rev. E 
61, 5278 (2000). 

J. Ophir and K. J. Parker, Ultrasound Med. Biol. 15, 319 (1989). 
^^J.-M. Escoffre, A. Novell, A. Zeghimi, and A. Bouakaz, Drug 
developpement and delivery 11 , 10 (2011). 

^^A. Novell, J.-M. Escoffre, and A. Bouakaz, Current Molecular 
Imaging 2, 77 (2013). 

^"^C. C. Church, J. Acoust. Soc. Am. 97, 1510 (1995). 

Hoff, P. C. Sontum, and J. M. Hovem, J. Acoust. Soc. Am. 
107, 2272 (2000). 

^®D. Chatterjee and K. Sarkar, Ultrasound Med. Biol. 29, 1749 
(2003). 

^^P. Marmottant, S. van der Meer, M. Emmer, M. Versluis, 
N. de Jong, S. Hilgenfeldt, and D. Lohse, J. Acoust. Soc. Am. 
118, 3499 (2005). 

^®M. Versluis, D. E. Goertz, P. Palanchon, 1. L. Heitman, S. M. 
van der Meer, B. Dollet, N. De Jong, and D. Lohse, Phys. Rev. 
E 82, 026321 (2010). 

Leighton, The acoustic bubble (Academic Press, London, 
2012) p. 640. 

20L. Ruby, Am. J. Phys. 64, 39 (1996). 

M. Braun and Y. S. Kivshar, The Frenkel-Kontorova model: 
concepts, methods, and applications, edited by Springer-Verlag 
(Springer, 2004) p. 472. 

^^B. Denardo, B. Galvin, A. Greenfield, A. Larraza, S. Putterman, 
and W. Wright, Phys. Rev. Lett. 68, 1730 (1992). 

^^W.-Z. Chen, Phys. Rev. B 49, 15063 (1994). 

^"^R. B. Thakur, L. English, and A. Sievers, J. Phys. D: Appl. 
Phys. 41, 015503 (2008). 

^^D. Armbruster, M. George, and 1. Oprea, Chaos 11 , 52 (2001). 
J. Sanchez-Morcillo, N. Jimenez, N. Gonzalez, S. Dos Santos, 
A. Bouakaz, and J. Chaline, Modelling in Science Education and 
Learning 6, 75 (2013). 

Sanchez-Morcillo, N. Jimenez, J. Chaline, A. Bouakaz, and 
S. Dos Santos, in Localized Excitations in Nonlinear Complex 
Systems (Springer, 2014) pp. 251-262. 




